The coherent scattering function of the reptation model: 
simulations compared to theory. 
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We present results of Monte Carlo simulations measuring the coherent structure function of a 
chain moving through an ordered lattice of fixed topological obstacles. Our computer experiments 
use chains up to 320 beads and cover a large range of wave vectors and a time range exceeding 
the reptation time. For additional information we also measured the coherent structure function of 
internal pieces of the chain. 

We compare our results (i) to the predictions of the primitive chain model, (ii) to an approximate 
form resulting from Rouse motion in a coiled tube, and (Hi) to our recent evaluation of the full 
reptation model, (i) The primitive chain model can fit the data for times t > 2OT2, where T2 
is the Rouse time of the chain. Besides some phenomenological amplitude factor this fit involves 
the reptation time T3 as a second fit parameter. For the chain lengths measured, the asymptotic 
behavior T3 ~ iV 3 is not attained, (ii) The model of Rouse motion in a tube, which we have 
criticized before on theoretical grounds, is shown to fail also on the purely phenomenological level. 
(Hi) Our evaluation of the full reptation model yields an excellent fit to the data for both total 
chains and internal pieces and for all wave vectors and all times, provided specific micro-structure 
effects of the MC-dynamics are negligible. Such micro-structure effects show up for wave vectors of 
the order of the inverse segment size and enforce the introduction of some phenomenological, wave 
vector dependent prefactor. For the dynamics of the total chain our data analysis based on the full 
reptation model shows the importance of tube length fluctuations. Universal (Rouse-type) internal 
relaxation, however, is unimportant. It can be observed only in the form of the diffusive motion of 
a short central subchain in the tube. 

Finally we present a fit formula which in a large range of wave vectors and chain lengths reproduces 
the numerical results of our theory for the scattering from the total chain. 



I. INTRODUCTION 



Dynamical properties of dense polymer systems like 
melts or dense solutions often are analyzed within the 
framework of the reptation model [1,2]. Reptation is a 
specific mechanism for the motion of a single tagged chain 
through an environment of other chains. It is based on 
the idea that the background chains act as impenetrable 
obstacles which confine the motion of the tagged chain to 
a tube roughly defined by its present configuration. The 
local motion of the inner parts of the chain is restricted 
to the diffusion of little wiggles of 'spared length' along 
the tube. Globally the motion is driven by the chain 
ends, where wiggles are created or destroyed. Creation 
of a wiggle shortens the tube by its spared length, de- 
struction prolongs the tube in some randomly choosen 
direction. In the long run this motion of the chain ends 
leads to the complete destruction of the original tube and 
to large scale diffusion of the chain. 

Formulated in more precise terms, the reptation model 
deals with the stochastic motion of a flexible chain em- 
bedded in a fixed environment of obstacles which form 
the edges of a regular lattice in three dimensional space. 
In this work we present results for the coherent struc- 
ture function measured in an extensive simulation of this 



model. The measured coherent structure functions of 
the total chain and of internal subchains are compared 
to the results of our recent analytical evaluation [3] of 
the model. For the total chain there exist previous ap- 
proximate theories based on reptation [4,5], which are 
included in the comparison. These theories do not treat 
the full dynamics of the model, but neglect so called 'tube 
length fluctuations'. For typical chain lengths used in 
(computer- or physical-) experiments these fluctuations 
are known to yield important contributions, as has first 
been pointed out in Ref. [6] in the context of an analysis 
of the viscosity. 

As mentioned above, results of the reptation model 
generally are used to analyze data for systems like poly- 
mer melts [7], where the surrounding of a given chain 
certainly is far from forming an ordered time indepen- 
dent lattice of obstacles. Clearly the surrounding chains 
slowly move away, which leads to 'constraint release' [8], 
an effect that becomes important [7] outside the limit of 
asymptotically long chains. Also disorder in the distri- 
bution of obstacles might lead to fluctuations in the local 
tube diameter, thus affecting the local mobility of spared 
length. In this work we omit all such effects of the en- 
vironment and study the coherent structure function of 
the pure reptation model, as described above. This is a 
necessary prerequisite for an analysis aiming at the sep- 
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aration of the different effects present in a real melt. 

To illustrate the problem we now briefly recall some 
typical results [f ,2] of reptation, as established for very 
long chains. We concentrate on the motion of an inter- 
nal segment, which from a theoretical point of view is the 
simplest quantity to discuss. 

Simple as it is, the reptation scenario involves several 
time scales and leads to a rich phenomenology, ft needs 
a microscopic time Tq before the chain feels the existence 
of constraints due to its surrounding. Generally T is 
taken as the Rouse time of a short subchain of N e seg- 
ments: To <~ N 2 , where the 'entanglement length' iV e is 
choosen such that the coil diameter of the subchain is of 
the order of the diameter of the tube, which substitutes 
the surrounding chains. The second time scale T 2 is the 
relaxation time of the total chain in a fixed tube, i.e. the 
time a wiggle needs to diffuse over the whole chain. It 
depends on chain length N as T 2 - T (N/N e ) 2 ~ N 2 
and thus behaves as the relaxation time of a free chain 
in the Rouse model. The longest scale T3 is the 'repta- 
tion time'. It measures the time which the motion of the 
chain ends needs to completely destroy the initial tube. 
In the limit of long chains the reptation model predicts 
[I] T 3 ~ (N/N e fT a . 

For observables like the motion of individual segments 
the model yields asymptotic power laws, where the ex- 
ponent depends on the time range. We here quote the 
results for the motion of the central segment j = N/2: 

I 5 \ f i 1/4 ; T Q « t « T 2 

( (r N/2 (t) - r N/2 (0)) ) ~ \ (t/N) 1 / 2 ; T 2 « t « T 3 
X 1 { t/N 2 ■ T 3 < t 

(l.f) 

We use the bar to denote the dynamic average, i.e. the 
average over the motion of spared length. The pointed 
brackets stand for the average over all initial configura- 
tions. 

Considerable effort has been invested to check these 
predictions in simulations of melts, but the outcome to 
date is not conclusive [7,9,10]. The t^-regime has never 
been properly identified. (Note that crossover from a 
region ((rj(t) — i"j(0)) 2 ) <~ t a , a < i, to free diffusion 
easily can pretend the existence of a t 1 / 2 -regime. What 
has to be demonstrated is the stability of this regime for 
a larger range of time and chain lengths.) Slowing down 
of segment motion in the range T < t < T 2 is observed 
[9], with an effective exponent somewhere between 1/4 
and 1/2. Only for some related observable, measuring 
the motion of a central segment relative to the center of 
mass, a i^-behavior seems to be established [10]. Real 
experiments do not measure {(vjit) — r, (0)) 2 ). However, 
a related quantity, the return-to-origin probability of a 
segment averaged over all segments, is measured in NMR- 
experiments. Here the equivalent to i^-behavior has 
been found in Ref. [11], but Ref. [12] reports equivalent 
results only for motion through a crosslinkcd gel, where 



constraint release is suppressed. The corresponding melt 
shows quite different behavior. 

Invoking tube length fluctuations and constraint re- 
lease we qualitatively may interpret the observed devia- 
tions from the asymptotic reptation results as crossover 
behavior outside the region of asymptotic chain lengths. 
However, there exist other theories of melt dynamics, 
which arc not based on the tube concept and describe 
many experiments as well [13,14]. (See also the review 
[7].) It thus is conceivable that the basic assumptions of 
the reptation scenario do not hold. To get more insight 
into these problems, we clearly have to quantitatively 
evaluate the consequences of the pure reptation model, 
beyond asymptotic power laws. 

In previous analytical work [15,16] we determined the 
motion of individual segments of the chain. Since our 
theory involves some approximations, we compared the 
results to simulations [17] of the Evans-Edwards model 
[18], which is an accurate implementation of the pure 
reptation model. In essence, both theory and simula- 
tions agreed in showing that the crossover among var- 
ious asymptotic power law regions is very slow. The 
crossover regions are so broad that the asymptotic power 
laws can be identified only for very long chains. For ex- 
ample, using the Evans-Edwards model with the small- 
est possible tube diameter we could identify the i 1 / 4 - 
law for the motion of the central segment only for chain 
lengths N > 160. This law is the easiest to observe, and 
our evaluation of the theory predicts that other asymp- 
totic laws unambiguously can be identified only for much 
longer chains. This result is in line with the known slow 
crossover behavior of the reptation time [6] , which is pre- 
dicted to reach the asymptotic law T 3 <~ N 3 only for 
chain lengths far beyond present day experimental fea- 
sibilities. Still, for the motion of individual segments, 
the full crossover functions can be calculated, and our 
analytical results very well agree with our simulations. 
Furthermore, for shorter chains all our analytical and 
simulational results qualitatively are very similar to re- 
sults of simulations of melts [9,10]. In later work [19] 
we considered chain motion in a time independent, but 
disordered environment, where the disorder affects only 
the chain mobility but leaves the equilibrium distribution 
of chain configurations unchanged. We found that with 
such 'kinematic' disorder reptation prevails, an observa- 
tion which recently has been supported [20] by rigorous 
bounds on the diffusion coefficient. Since such kinematic 
disorder certainly is present in a melt, these results again 
support reptation as adequate theory of melt dynamics. 

In contrast to the motion of individual segments the 
coherent structure function S c (q,t;N), (q = |q| : scat- 
tering vector, t : time), is well accessible in real experi- 
ments. It can be measured, for instance, by neutron scat- 
tering from a mixture of deuterated and hydrogenated 
chains. As mentioned above, approximate asymptotic 
forms for S c (q,t; N), based on reptation type theories, 
can be found in the literature [4,5]. However, as for the 
motion of individual segments we can evaluate the full 



2 



reptation theory for S c (q,t;N) also outside the asymp- 
totic regime. As for individual segments, we then expect 
to see important preasymptotic or crossover effects. The 
evaluation of S c (q, t; N) including full reptational dynam- 
ics is somewhat complicated, and our theory in detail has 
been presented in Ref. [3]. Here we compare the results 
to simulations of the Evans-Edwards model. We consider 
both the scattering from the total chain and from inte- 
rior subchains. The latter is important to estimate the 
reliability of the theory, which can treat end-effects only 
in some approximation. 

In Section 2 we discuss our simulations. Section 3 is de- 
voted to a comparison with the expressions for S c (q, t; N) 
given by Doi and Edwards [4,2] or by de Gennes [5], re- 
spectively. In Section 4 we qualitatively describe our the- 
ory and give an empirical fit formula which reasonably 
well describes our quantitative results. The comparison 
between our theory and our Monte Carlo data is pre- 
sented in Sect. 5. Section 6 contains our conclusions. 



(The imaginary part 5(q, i)C(q, 0) - C(q, i)«S(q, 0) of S c 
is zero on average, of course.) To get more information 
on the internal motion, we also measured the coherent 
structure function S c (q, t^ MC ^ ; M, N) of the central sub- 
chain of length M, defined by restricting the summa- 
tions in Eqs. (2.1), (2.2) to the interval [(N - M)/2 + 
1,(N + M)/2]. ^From our previous work, we expect 
to see features characteristic of reptation for N > 100, 
and we therefore used chain lengths N = 80, 160, 320. 
Since for our model, the entanglement length is esti- 
mated as N e 3.7 (sec Sect. 4.B), this yields values 
22 < N/N e < 87. Similar values are extracted from 
many simulations or experiments on melts, so that our 
results should be relevant also for the interpretation of 
such data. Monte Carlo time f( MC ) is measured in units 
of one attempted move per bead on average. We per- 
formed runs up to t^[^ = 5 • 10 10 , and wc measured the 
structure function up to i( MC ) = 4.5 • 10 9 using a moving 
time average. 



II. SIMULATIONS 



B. Statistical accuracy 



A. The Evans-Edwards model 

A very simple model for simulating reptation has been 
introduced by Evans and Edwards [18]. The configura- 
tion . . . , rjv} of the chain is taken as a random walk 
of N— 1 steps \rj — Tj-i\ — £o on a cubic lattice. The lat- 
tice spacing £ henceforth defines the unit of length. The 
obstacles are taken as the edges of the dual lattice. In the 
interior of the chain, the obstacles suppress any motion 
except for the motion of 'hairpins', i.e., configurations of 
three subsequent beads of type {rj-i, rj, tj + i = rj_i}. 
In an elementary move the tip rj of the hairpin with equal 
probability hops to any of the six neighbors of the site 
rj-i = tj + i. The chain ends ri,rjv are free to hop be- 
tween all neighbors of T2, rjv-i, respectively. Fig. 1 shows 
a sequence of internal configurations resulting from this 
dynamics. In our simulations we used the same imple- 
mentation of the model as in our previous work [17], and 
we measured the coherent structure function defined as 



N 



where 



= (C(q, t(^))C(q, 0) + 5(q, t( M ^))5(q, 0) 

(2.1 



N 



C(q,t) = 2cos(qr j (t)) 



N 



<S(q,t) =X) sin ( qr ^)) • 



(2.2) 



During a run, the longer chains do not diffuse very 
far, and data of a single run therefore are strongly cor- 
related. To get reasonably accurate results we have to 
average over many independent runs. A priori, this poses 
a problem for larger momenta, qR g > 1 , where R g is the 
radius of gyration of the chain. It is easily checked that 
the reduced width of the distribution of the static struc- 
ture function with increasing qR g rapidly tends to 1 . In- 
deed, in the limit N — > oo, q > fixed, the distribution of 
S c (q, 0; N) approaches the exponential distribution. This 
suggests that 10 4 runs are needed to reduce the statis- 
tical uncertainty to a few percent. This would pose no 
problems, if we just were interested in static properties. 
However, for the longer chains a single run extending to 
times well beyond the reptation time takes several hours 
on a standard work station. In measuring dynamic quan- 
tities the number of runs therefore inevitably is much 
smaller than needed for a precise determination of static 
quantities. 

Fortunately it turns out that the dynamics es- 
sentially is decoupled from the static configuration. 
This is illustrated in Fig. 2, which shows results for 
S c {q,t( MC 1;N),qlo = 0.5, TV = 320, normalized to the 
exact static structure function So(q,N) of the model, 
x which easily is calculated analytically (see Sect. 2.D). 
Each curve in Fig. 2 a is averaged over 10 3 independent 
short runs 

( t (MC) < 1Q 5) 5 inc i uding the 

moving time av- 
erage for each run. Clearly the scatter of S c (q, 0; N) 
is consistent with the above discussion. It is larger 
than the temporal variation of the curves. However, 
plotting in Fig. 2 b the normalized time dependence 
[S c (q, t^ MC ^\N) - S c (q,0;N)}/S o (q,N), we find that all 
curves nearly coincide. The global dynamics measured by 
the structure function, is not correlated with the static 
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configuration, an observation which supports one of the 
basic assumptions of the reptation model. The reason 
behind that observation is easily understood. The mo- 
bility of the chain is governed by the number of hairpins 
which essentially is Gaussian distributed and fluctuates 
much less than the static structure function. Further- 
more, except for rare events, viz. extremely stretched or 
extremely compact configurations, the number of hair- 
pins is independent of the overall (tube-) conformation 
of the chain. With this insight, we take as basic data 
the difference S c {q,t^ MC ^;N) - S c (q,0;N) for each run. 
The error bars in our plots give twice the standard devi- 
ation of this difference. We always plot the normalized 
dynamic structure function defined as 



S c (q,t 



(MC). 



N) = 1 



1 



(2.3) 



S {q,N) 

(s c (q,t^;N)-S c (q,0;N)) 

V V / / averaged over runs 

where So(q, N) is the exact static structure function, not 
the measured average value of S c (q,0;N). Depending 
on the number of independent runs, these two quantities 
differ by 0.1 - 6%. 

For reasons of computer memory, we performed runs 
with three different values of the maximal time twx\ 
To observe the small initial effects we for each value of 
chain length N and wave vector |q| averaged over 10 4 

short runs with t^x^ = 10 5 , taking data up to £( MC ) = 
1.4- 10 4 . These data are denoted by small dots in the later 



figures. In the intermediate time range, i£Jf£^ = 10 8 



measured the structure function for 10 4 < £( MC ) < 10 7 
and always performed 50 runs (big dots in the figures). 
In the long time range, — 5 • 10 10 , we took data 

for 10 5 < i( MC ' < 4.5 • 10 9 , averaging over 30 to 100 
runs (circles in the figures). For each set of runs with 
given maximal time, we also measured the average of 
S c (q,Q;N). 

For additional information on the statistical accu- 
racy of our data, we measured the imaginary part of 
S c (q, t( MC ^ 7 N), which rigorously vanishes for t^ MC ^ = 0, 
but fluctuates about zero for £( MC ) > in a finite sam- 
ple. For longer times we typically found average values of 
order ±0.01So(q, N), again smaller than the uncertainty 
of S c (q, 0; N) in long-time runs. 



C. Momentum range 

First considerations suggest to restrict the analysis to 
the momentum range R~ 2 « ? 2 « £q 2 . Momenta of 
order q 2 R 2 g ^ 1 do not resolve the tube but rather see 
a cloud of beads. For momenta q 2 £\ > 1 the micro- 
structure of the chain might play a role. In our simu- 
lations the above condition can only poorly be satisfied. 
For our longest chain (N = 320) it reads 0.02 < q 2 i\ < 
1, leaving at best a small window close to qlo ~ 0.4. 



However, a closer inspection reveals that these con- 
siderations do not seriously restrict dynamic measure- 
ments. To check the relevance of the condition q 2 l\ <C 1 
we analyzed the dynamics of a free chain. As is well 
known, asymptotically the standard Rouse dynamics is 
found if for a lattice chain 'kink jumps': {ij — rj_i = 

Sl, r j+ i - Y 3 = S 2 } {Vj - Vj-! = S 2 , Ij + i - Vj = 

Si}, Si • s 2 = are allowed besides hairpin- moves. We 
found that with this modified dynamics, the normalized 
coherent structure function is in excellent agreement with 
the relaxation function calculated from the continuous 
chain Rouse model, even for q£o = 1. Micro-structure 
effects arising from the finite segment size die out very 
rapidly on the scale of about 10 MC-steps and there- 
fore should be negligible also for reptation dynamics. Of 
course this does not exclude the possibility of dynamical 
micro-structure effects which might influence the short 
time regime and are not accounted for by the reptation 
model. In particular, reptation docs not properly treat 
the dynamics of fluctuations perpendicular to the tube 
axis. 

Concerning the condition q 2 R 2 ^> 1, we note that 
smaller wave vectors, of course, provide little informa- 
tion on the details of the internal motion of the chain, 
but at least they measure the global motion of the coil. 
However, the slowness of this motion asks for extremely 
large time ranges t^£P 3> T3. We therefore performed 
only one series of simulations for |q|^o = 0.1, N = 160, 
corresponding to q 2 R 2 g = 0.267. Most of our simula- 
tions use values 3 < Q 2 R^ ^5 50, with a maximal value of 
.202 



we q 2 R 2 w 53 reached for \q\£ = 1, N = 320. 



D. Normalization 

In our simulations we choose q to point into one of 
the three lattice directions. For that choice, the static 
correlations between segments j and k take the form 



o iq(r :j (0)-r Ji (0))\ = ( £ 1 



where 



= - In 



3 ' 3 C ° S |q| 
cxp(-q 2 \k - j\) 



cos q 



\k-j\ 



(2.4) 



(2.5) 



Recall that the lattice constant £q defines the unit of 
length. Summing the segment indices over the chain, we 
find the static structure function which is used in nor- 
malizing our results: 



S (q,N) = CX P(-'/ 

j;k = l 

= N + 2 



(2.6) 



(e<? 2 - l) 2 



e -f{N-i) _ N +{N _ l)e f 
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In our simulations, we for each run averaged over all three 
lattice directions. 



III. COMPARISON OF DATA AND SIMPLIFIED 
REPTATION TYPE THEORIES 

A. Asymptotic form of S c {q,t;N) derived by Doi and 
Edwards 

A simplified version of the reptation model concen- 
trates on the dynamics of the 'primitive chain' [2,4], 
which is a reduced form of the chain, lying stretched 
in the tube. In the Evans-Edwards model the primitive 
chain can be viewed as the non-reversal random walk de- 
rived from the random walk configuration of the physical 
chain by cutting off all hairpins. All internal degrees of 
freedom are neglected so that within time interval At, 
all parts of the primitive chain move the same distance 
As along the tube. The length of the primitive chain is 
taken to be fixed. (This model often is addressed as 'the 
reptation model'. We will use the term 'primitive chain 
model' to distinguish it from the full reptation model 
which deals with the dynamics of spared length as the 
elementary process.) 

With its simplifications, the primitive chain model 
treats only the destruction of the initial tube, as result- 
ing from the global motion of the chain. It neglects tube 
length fluctuations which are due to the uncorrelated mo- 
tion of the chain ends as well as relaxation in the inte- 
rior of the tube. Since both these additional effects are 
governed by the equilibration time T 2 of the chain, the 
primitive chain model is restricted to times t 3> T 2 . Fur- 
thermore, the segment indices are taken as continuous 
variables, and this 'continuous chain limit' restricts the 
theory to long chains: N ^> 1. 

Within this model the time dependent correlation func- 
tion of two beads: 



and 



S(q,t;j,k,N) = (exp[iq(r,(t)-r fc (0))] 



(3.1) 



obeys a diffusion equation. Solving this equation and in- 
tegrating j, k over the chain, Doi and Edwards arrive at 
the following result for the normalized coherent structure 
function (see Ref. [2], chapter 6.3.): 



S c (q,t;N) 



Sde(Q,t) 



Sc(g,t;N) 
S c (q,0;N) 



S DE (q 2 R 2 g ,t/T 3 ) (3.2) 



Q sin 2 a p 



D(Q)^a 2 JQ 2 /4 + Q/2 + a 2 p ) 



exp (-apr) 



(3.3) 



Here the a p are the positive solutions of 



D(Q) = -^(e- Q -l + Q) 



(3.5) 



a p tan a p = \Q= \l 2R l ' 



(3.4) 



is the Debye function. With the assumptions of the the- 
ory, the reptation time attains its asymptotic behavior 
T 3 <~ TV 3 . Note that T 3 is related to the time scale r<j 
introduced in Ref. [2] by T 3 = \t&. 

The result (3. 2), (3. 3) for S c (q,t;N) depends only on 
Q = q 2 R 2 g ~ q 2 N and r = t/T 3 ~ i/TV 3 , but not on TV 
separately. To test this feature, we carried through sim- 
ulations for TV = 320, q = 0.25 and TV = 80, q = 0.5024, 
both parameter sets resulting in q 2 R 2 ~ 3.3. Fig. 3 com- 
pares our simulation results to Sde, plotted as function 
of log 10 (t/T 3 ). The Monte Carlo time has been scaled so 
that in the region t/T 3 > 1, the data for TV = 320 fit to 
the theoretical curve. In view of T 3 <~ iV 3 , for TV — 80 an 
additional factor 4 3 has been included in the time scale. 
With this scaling, the deviation between the two sets of 
data shown in Fig. 3 proves that we have not yet reached 
chain lengths large enough for the primitive chain model 
to hold. 

Still, for larger times, Sde and the data for TV = 320 
agree quite well. This suggests to treat T 3 = T 3 (N) as a 
fit parameter in adjusting the data to SDE{q 2 R%t/T 3 ), 
giving up the strict proportionality T3 <~ TV 3 . However, 
a closer inspection of Fig. 3 reveals that the data ini- 
tially decrease faster than Sde- This is a systematic 
effect, observed for all chain lengths and wave vectors. 
It points to the influence of relaxation modes neglected 
in the primitive chain model. According to a suggestion 
of de Gennes [5], internal relaxation, in particular, yields 
an initial decrease of S c (q, t; TV) which saturates at some 
q-dependent plateau value. We thus should fit the data 
to BDESDE{q 2 R 2 g ,t/T 3 ) 1 with Bfl£ as another free pa- 
rameter. 

A detailed analysis of the reptation model points to 
tube length fluctuations as the origin of the initial de- 
crease, rather than internal relaxation. (See Sect. 4.C 
and Ref. [3].) Still we may fit our data to BdeSde, 
where Bde = B DE (q,N), T 3 = T 3 (TV). Fig. 4 shows 
the results for TV = 320. We clearly find very good 
agreement among theory and data in the time region 
t > 20T 2 . (The estimate for T 2 has been taken from 
our theory, see Sect. 4.B.) Deviations occur for smaller 
times, which is consistent with the approximations in- 
herent in the primitive chain model. Within the realm 
of that model the central segment moves according to 
• ((r N /2(t) - r N/2 (0)) 2 ) ~ (t/TV) 1/2 ,T 2 « t « T 3 (cf. Eq. 

(1.1)), and from our previous work [17] we know that this 
law certainly is not attained before t > 20X 2 . The results 
shown in Fig. 4 are typical also for other chain lengths. 

Since T 3 now plays the role of an effective fit parameter, 
it a priori could depend both on TV and q. Parameters 
T 3 = T 3 (q,N), B DE = B DE (q,N) extracted by fitting 
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our data are collected in Table 1. Any g-dependence of 
T 3 is found to be weak, if significant at all. Boe depends 
on q, but is essentially independent of N. 

Ignoring any g-dependence, we in Fig. 5 have plotted 
the values of T 3 (iV) normalized to the asymptotic be- 
havior T 3 °° resulting from Eqs. (4.5), (4.7), (4.12) be- wher e 
low. It must be stressed that we here define a reptation 
time in terms of the large time behavior of the scatter- 
ing function (and the Doi-Edwards form Sjje)- Other 
definitions based on other observables (or other theoreti- 
cal expressions) may yield somewhat different results. In 
previous work [17], we defined a reptation time T 3 by 
the criterion that the mean squared distance moved by 
a chain end equals the equilibrium mean squared end- 
to-end distance: ((r (T 3 ) -r (0)) 2 \ = R\. The previ- 



model of a Rouse chain in a coiled tube yields a 'local' 
contribution 



ous data are included in Fig. 5 (open circles). The two 
definitions yield somewhat different results for the cor- 
rections to the asymptotic behavior. In the range of N 
measured here, T 3 is about 30 % below T 3 as taken from 
the scattering function. It, however, must be noted that 
the values extracted from the scattering functions depend 
somewhat on the time range included in the fit. We es- 
timate this uncertainty to be of the order of 5 %. In 
Fig. 5 we included lines corresponding to effective power 
laws T 3 ~ N z . The effective powers are consistent with 
expectations based on previous work [6] . 

To summarize, our results show that the coherent 
structure function of the primitive chain model may well 
be used to fit data for large times, t > 20T 2 , provided we 
allow for some phenomenological prefactor Bde and take 
the reptation time T 3 as an effective parameter defined 
by the fit. For smaller times, deviations are seen, that 
increase with increasing q. 



B. Comparison to Gennes' theory 

In his work [5], de Gennes considers only intermediate 
values of q: £q ^> q 2 3> Rg 2 and constructs the scat- 
tering function as a sum of two terms. For t » T 2 the 
'creep' term dominates. This term is just the limiting 
form of S DE (Eq. (3.3)) for Q -> 00: 



5W(f,JV)4^(2p+l)- 2 exp 



p=0 



IT 2 t 



-<* + 1 >T5i 



(3.6) 



It tends to 1 for t/T$ — ► 0. It is combined with some 
contribution of local relaxation, which is calculated in 
two steps. First the interior relaxation of a Rouse chain 
in one dimensional space is calculated, where the chain 
is stretched so that the end-to-end distance equals the 
tube length. In the second step this one dimensional 
motion is embedded into the three dimensional random 
walk configuration of the tube. Tube length fluctuations 
are neglected. With some additional approximation, this 



S w (q,t) = e* 1 erfc V*T 



(3.7) 



(3.8) 



Since again the validity of the asymptotic law T 3 <~ 
N 3 /N e is assumed, the chain length N and the entan- 
glement length N e drop out in Eq. (3.8), as expected for 
a contribution resulting from strict one dimensional in- 
ternal relaxation. (q, t) describes internal relaxation 
on length and time scales exclusively determined by q. 
The total result for the normalized coherent scattering 
function reads 



S c (q,t;N) = S dG (q,t;N) = 

(l-B dG )S^(q,t) + B dG S^\t,N) 



(3.9) 



where according to de Gennes 



B dG = B dG (q) = 1 



^q 2 £ 2 

36 q 







(3.10) 



is independent of N. 

In an attempt to extend the range of wave vec- 
tors to q£o > 1, Schleger et al [21] used B dG (q) = 
exp(-_q 2 £lN e /36). With this choice it was found [21-23] 
that S dG (q, t; N) over a large range of wave vectors and 
for several chain lengths yields a good fit to neutron scat- 
tering data from melts of Polyethylene. However, Piitz 
et al [10], using the same form of S dG (q,t; N) to ana- 
lyze simulation data for melts of very long chains, found 
only poor agreement. In later work [24] they argued that 
this was due to some ambiguity in the relation between 
entanglement length N e and tube diameter. 

In our recent work [3] we reconsidered the model of a 
Rouse chain in a tube. Our analysis reveals some serious 
deficiency of this model: it does not start from equilib- 
rium initial conditions. An ensemble of stretched one 
dimensional random walk chains folded into the three 
dimensional random walk configuration of the tube is 
not identical to the equilibrium ensemble of three di- 
mensional random walk chains. The local structure of 
the chains differs on the scale of the tube diameter. 
For the static structure function, this yields a correc- 
tion of relative order N e /N that vanishes in the limit 
q 2 R 2 g = const, N — ► 00. For the time dependence, how- 
ever, the effect is serious since the non-equilibrium initial 
conditions relax only on scale T 2 . Our analysis shows that 
this unphysical relaxation indeed dominates the time de- 
pendence of the 'local' contribution (q, t) , as calcu- 
lated from this model (see Fig. 3 of Ref. [3]). Data analy- 
sis based on Eq. (3.9) with (q, t) taken from the model 
of a Rouse chain in a tube therefore is not particularly 
meaningful. 
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Still, in view of the experimental findings cited above, 
we may ask whether Eqs. (3.6) - (3.9) can be used as 
heuristic modeling of the coherent structure function 
resulting from reptation. Fig. 6 shows the result for 
the longest chain (N = 320) and largest wave-vector 
(q = 1.0) measured. The resulting value q 2 R 2 g ~ 53.3 
is large enough for a reasonable test of the form (3.9), 
which assumes q 2 R 2 g 3> 1- Replacing in Eq. (3.9) S c by 
the full result Sde of the primitive chain model does 
not seriously change the picture. The reptation time is 
fixed by the long-time tail and is taken from Table 1. 
For BdG, we used the form suggested by Schleger et al: 
BdG = cxp(— q 2 if,N e /36). N e in principle is known from 
our previous analysis of segment motion: N e = 3.69 (see 
Sect. 4.B). Thus all parameters are fixed, and the re- 
sult for SdG (full curve in Fig. 6) strongly deviates from 
the data, except for the extreme long-time tail. This 
is no surprise since the value BdG ~ 0.9, resulting for 
N e = 3.69, considerably exceeds the value Bde ~ 0.78 
extracted from fitting with the Doi-Edwards form. (See 
Table 1.) If we treat BdG as a free parameter, we clearly 
in the range t > 2OT2 can enforce agreement among the- 
ory and data, (see Fig. 4), at the expense of considerably 
underestimating S c for t < Ti. The situation can be im- 
proved only, if we also change the scale of t\, dividing t\, 
(Eq. 3.8), by a factor of order 200. Quite similar results 
are found for N = 160, where t\ has to be divided by a 
factor of order 50 to reproduce the average trend of the 
data for t < T 2 . 

We conclude that also from a purely phcnomcnologi- 
cal point of view, the form (3.9) of S c is not justified. 
The large and chain length dependent rescaling required 
for t\ suggests that tube length fluctuations governed by 
time scale Ti might be much more important than inter- 
nal relaxation processes governed by an iV-independcnt 
scale. 

We finally note that in Ref. [23] a modified form of Eq. 
(3.9) has been proposed, in which the creep term is 
calculated from tube length fluctuations, but the form 
(3.7) of the local contribution is retained. The analy- 
sis is restricted to the time range t -C Ti, which is the 
relevant range for neutron scattering experiments. From 
Fig. 6 it is clear that with the proper parameter values 
N e = 3.69, t\ as given in Eq. (3.8), this modification 
cannot improve the fit for t <^ Ti, since tube length fluc- 
tuations only decrease the contribution S^ c \ Again a 
reasonable fit can be reached only if we take both BdG 
and the scale of t\ as free parameters. 



IV. RESULTS OF THE FULL REPTATION 
MODEL 

A. Basic ideas of our approach 

The primitive chain model neglects all internal degrees 
of freedom, and an attempt to model internal relaxation 



as one dimensional Rouse motion leads to unphysical ini- 
tial conditions, which seriously affect the results up to 
times of the order of the internal relaxation time Ti. 
However, the simplifying assumptions underlying these 
approaches are no essential part of the original reptation 
model. As recalled in the introduction, reptation as sin- 
gle elementary process involves the diffusion of spared 
length along the chain, together with its decay and cre- 
ation at the chain ends. The separation of the dynam- 
ics into internal relaxation, tube length fluctuations and 
motion of the primitive chain therefore is somewhat ar- 
tificial. In particular, tube length fluctuations cannot 
properly be separated from the motion of the primitive 
chain, as will be discussed in Sect. 4.C. Furthermore, 
from our analysis of segment motion we know that for 
typical experimental chain lengths we are in a crossover 
region where all dynamic effects must be treated on the 
same level by evaluating the consequences of full repta- 
tional dynamics. 

In our analytical work [16,3] we use a very simple im- 
plementation of the reptation idea, in which the wiggles 
of spared length are represented by particles which hop 
along the chain, with hopping probability p per time step. 
The particles do not interact, and a given particle sees 
the others just as part of the chain. If a particle passes 
a bead, it shifts the position of this bead in the tube 
by the spared length is- The chain ends are coupled to 
large particle reservoirs, which absorb and emit particles 
at such a rate that the equilibrium density po of spared 
length on the chain is maintained on the average. Ab- 
sorption or emission of a particle prolongs or shortens 
the tube at the corresponding chain end by the spared 
length is- Keeping track of the change in the occupation 
number of the reservoirs, we therefore control the tube 
length fluctuations as well as the destruction of the orig- 
inal tube. In particular, within time interval [0, t] an end 
piece of length isn max (t) of the original tube has been de- 
stroyed, where n max (t) is the largest negative fluctuation 
of the occupation number of the corresponding reservoir 
during this time interval. 

All moments of the stochastic processes which deter- 
mine the motion of internal segments or the occupation 
of the reservoirs, can be evaluated rigorously, but the de- 
termination of the maximal fluctuation n max (t) poses a 
serious problem. The occupation number of a reservoir 
carries out a correlated random process, since a particle 
emitted can be reabsorbed by the same reservoir later. 
This correlation dies out only on scale of the internal 
equilibrium time Ti of the chain. It does not prevent the 
evaluation of arbitrary moments, but the maximal fluctu- 
ation cannot be calculated rigorously. Such a calculation 
is possible [25] only for an uncorrelated random process. 
To determine the degree of tube destruction, we therefore 
have to resort to some approximation. In our method, 
we basically for each final time t replace the correlated 
random process by that uncorrelated random walk which 
for this time yields the correct moments. The effective 
hopping rate of this random walk depends on the final 
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time t. It changes from the microscopic rate po ■ p for 
t <C Ti to the mobility po-p/N oi the primitive chain for 
t>T 2 , which is a physically most reasonable behavior. 

In essence, this 'mean hopping rate' approximation for 
the coherent structure function smoothly interpolates be- 
tween two rigorously accessible limits. For i < T2, tube 
renewal does not influence the motion of an interior piece 
of the chain, and the coherent structure function of such 
pieces can be calculated rigorously. For t ^> Ti the corre- 
lations of the stochastic processes are negligible and the 
mean hopping rate approximation should become exact. 
Indeed, we find that in this limit the motion of all chain 
segments is tightly bound to the motion of the chain ends. 



Here dtoV*(jo, to\m), m = 0, N, is the probability that 
the initial tube is finally destroyed within time interval 
[to, to + dto], its last piece being the initial position of 
segment j , which at time to is occupied by chain end 
or N, respectively. The exponential factors result from 
the random walk configuration of the tube, cf. Eqs. (2.4), 
(2.5). The inhomogeneity S tyT \q,t;j,k,N) is the contri- 
bution to S(q, t; j, k, N) of all those stochastic processes, 
which do not completely destroy the original tube. 

Summing Eq. (4.2) over the segments j, k we find a 
system of two equations which have to be solved numeri- 
cally. The kernel V* and the inhomogeneity are de- 
termined within the mean hopping rate approximation. 
The results are lengthy and will not be reproduced here. 
We rather in Subsect. D give an analytical expression, 
which in the range of wave vectors and chain lengths con- 
sidered in the present work, reasonably well reproduces 
the numerical results of our theory. 

We finally note that we analytically can prove (Ref. [3] , 
Sect. 7) that our theory in the limit of infinite chain 
length N — > co, with t/T$ and q 2 R 2 g kept fixed, repro- 
duces the result of the primitive chain model. Also the 
relation to the model of a Rouse chain in a tube can 
be analyzed in precise terms, if we consider an interior 
piece of length M in an infinitely long chain. We find 
(Ref. [3], Sect. 5) that this model reproduces the results 
of reptation only for t/T 2 (M) -> oo, q 2 R 2 g (M) fixed, 
where T 2 (M) or R g (M) are the equilibration time and 
the radius of gyration of the subchain considered. For 
tjT^iM) 1 the non-equilibrium initial condition seri- 
ously affect the scattering function, as has been discussed 
in Sect. 3.B. 



As a result, the problem reduces to the uncorrelated mo- 
tion of a single stochastic variable, as in the primitive 
chain model. The details of our approximation are dis- 
cussed extensively in Refs. [3,16], and will not be repeated 
here. 

The coherent scattering function can be determined by 
summing the contribution of two beads 

S(q, t; j, k, N) = (^^JW^AP))^ (4.1) 

over the bead indices j, k. For this function we in Ref. [3], 
Sect. 6. A, have derived an integral equation of the form 



(4.2) 



B. Microscopic parameters of the reptation model 

The microscopic parameters of the model are the seg- 
ment size £o = r j-i I , the average density po of mobile 
particles on the chain, the spared length per particle is, 
and the hopping rate p. Measuring all lengths in units of 
lo, we introduce the dimensionless spared length 

h = l f • (4.3) 

It turns out that for all times beyond truly microscopic 
times t < 2/p, the hopping rate p combines with t to yield 
the time variable 

i = pt , (4.4) 

which we will use in the sequel. The relation of t to the 
Monte Carlo time t^ MC ^ introduces the fit parameter To: 

t = T t^ . (4.5) 

Also po and is combine into a single important parame- 
ter. The number of particles that passed over a bead on 
average increases with time, and if it is sufficiently large, 
the discreteness of the individual hopping processes be- 
comes irrelevant. The results then depend only on the 
combination i\po- In our previous work on the motion 
of individual segments [15-17], we found that is and po 
separately enter the results only for t < 10 3 . In this 
time region, a segment on average has moved less than 
10 steps in the tube and still feels the discreteness of the 
process. For the structure function we find that not even 
this small time region is seriously affected. If we ignore 
the discreteness of the hopping process, the results for all 
t > 1 change by less than 0.5 %, which coincides with the 
accuracy of our numerical evaluation. Thus the only mi- 
croscopic parameters relevant for the coherent structure 
function are the combination £ s po and the time scale tq. 



S(q,t;j,k,N) = S^ T \q,t;j,k,N) 

+ / dt Q J2 \V*(jo,to\0)exp(-q 2 \jo-k\)S(q,t-to;j,0,N) 
Jo io=o 1 

+ V*(jo,to\N)exp(-q 2 \jo -k\)S(q,t-t ;j,N,N)} . 
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The reason behind this empirically observed suppression 
of initial discreteness effects will be discussed in Subsect. 
4.C. 

We furthermore note that also the discreteness of the 
underlying chain turns out to be unimportant. Pro- 
vided we normalize the coherent structure function by the 
static structure function calculated for the same micro- 
structure, we for N > 50 and wave vectors qlo < 2 within 
the above quoted accuracy find the same results for a 
continuous chain as for the model where we sum over 
discrete bead indices j, k. This observation is consistent 
with the independence of statics and dynamics discussed 
in Sect. 2.B and illustrates that our calculation indeed 
yields universal results. 

Since in our simulations we use the same model as in 
our previous work on segment motion [17], we can take 
the numerical values of the microscopic parameters from 
there. In analyzing the Monte Carlo data, we thus use 
the value 



P sPo = 1.23 



(4.6) 



Having a much larger set of data available than previ- 
ously, we somewhat readjust the time scale. We use 



r = 6.8 • 10" 



(4.7) 



rather than the previous value t = 6.092 • f 0~ 2 . For the 
logarithmic scale log 10 t of the figures in Ref. [17], this 
amounts to a shift by —0.048, which does not change the 
good agreement among theory an experiment which in 
Ref. [17] is shown to hold over about 6 decades of time. 
In our formulation of the reptation model, the entan- 

1/2 

glement length N e , or the tube diameter £oN e , equi- 
valently, do not show up explicitly. They are hidden in 
the spared length and the density of the particles, i.e., in 
the overall mobility of the chain. In contrast, the previ- 
ous approaches explicitly involve these parameters. As 
discussed at the end of the last subsection, our theory in 
appropriate limits reproduces the previous results, which 
allows us [3] to relate our parameters is,po,P to those 
of these other models. Specifically, we find a relation for 
the entanglement length: 



N e = 3P sPo 



(4.8) 



With the numerical value (4.6) this yields N e = 3.69 
and implies that for the Evans-Edwards model with the 
smallest possible obstacle spacing it needs about 4 steps 
before the obstacles come into play seriously. 

We note that IsPo an d thus N e have been determined 
[17] by fitting data for the motion of individual segments 
in the time range t > 10 3 . N e is thus not influenced by 
the deviations from the full reptation model occuring in 
the initial range t < 10 3 , as discussed below, (Sect. 5). 
Thus consistency of the analysis clearly enforces use of 
the same value N e for all observables. We estimate the 
uncertainty of this value to be in the range of 5 %. 



We now also can give a quantitative definition of the 
time scales. Identifying the equilibration time T 2 with 
the Rouse time of a free chain of N segments we from 
the asymptotic relations among the models find 



T 2 = pT 2 = 



(N + lf 



(4.9) 



The reptation time can be defined as the average life time 
of the original tube 



POO 

f 3 = P dt t r*(t Q ) . 

Jo 



(4.10) 



Here dt P*(to) is the probability that the tube is finally 
destroyed within time interval [to , to + dt ] . It is related 
to the distribution V*(jo, to\m) introduced in Eq. (4.2) 
via 



N 



V*(t )= [V*(j ,t \0) + V*(j ,t \N)\ ■ (4.11) 

jo=0 



Evaluation in the asymptotic limit N — > oo yields 



iV 3 3 7V 3 



4P sPo AN e ■ 



(4.12) 



In this limit, f 3 is related to the reptation time Td intro- 
duced for the primitive chain model [2] as 



j-(oo) = ^_ 



PTd , 



(4.13) 



which is the relation used in Sect. 3 in our analysis of the 
primitive chain model. 



C. Qualitative discussion of the coherent structure 
function 

Tube destruction, tube length fluctuations, and inter- 
nal relaxation are different facets of the same microscopic 
dynamics. Therefore any separation of these processes 
is to some extent artificial. Yet some rough discussion 
based on these concepts is useful with regard to the in- 
terpretation of the data presented in Sect. 5. 

We first consider tube length fluctuations AL. Here 
a precise definition is possible by identifying AL with 
the fluctuations of the total spared length. On average 
the spared length equals IsPoN, where p N is the av- 
erage number of particles on the chain. Since the fluc- 
tuations of the particle number essentially are Gaussian 
distributed we find 



AL - , 

-j-nesVtoN 



-N e N 
3 



(4.14) 



where N e is introduced via Eq. (4.8). Within the equi- 
libration time T 2 tube length fluctuations on average 
on both chain ends replace a piece of length AL/2 of 
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the original tube by a new piece. In real three dimen- 
si onal spac e, this new piece has an average extension of 
\JIqKL~/2, and to observe the effect, the scattering vector 
must obey q 2 £ AL/2 > 1, or 



M 



n 2 f 2 > 



12 



(4.15) 



equivalcntly. For large wave vectors, the coherent struc- 
ture function S c (q, t; N) is determined by the part of the 
original tube that still exists at time t. Taking only tube 
length fluctuations into account, we therefore find for 



S c (q,T 2 ;N) = 



S c {q,t;N) N-AL/e 




(4.16) 



provided q is large compared to the bound (4.15). This 
estimate yields the effect of fully developed tube length 
fluctuations. With increasing chain lengths it clearly is 
suppressed, but it should be well visible up to fairly long 
chains N/N e w 10 2 . 

Tube destruction in the sense of the primitive chain 
model is responsible for the main part of the decay of 
the coherent structure function. According to Eq. (4.16), 
it is the dominant effect as soon as 1 — S c (q,t,N) 3> 
a/ N e /3N. For t —* T 2 , it is strongly correlated with 
tube length fluctuations which, seen from a microscopic 
point of view, drive the tube destruction. It therefore is 
no surprise that evaluating the coherent structure func- 
tion of the primitive chain model, §DE{q 2 R 2 g ,t/T 3 ), for 
q 2 R 2 g ^> 1 and time t = T 2 , we find a decrease of the 
same order of magnitude as that due to tube length fluc- 
tuations. 



I ( 'riff, - tjt 



l-S (c) (T 2 ,A0 



- Tl + (9^-const T 3 /T 2) _^_^e, Q{ -const Ns 

(4.17) 



Recall that S^(t,N) (Eq. (3.6)) is the limit of S DE for 
q 2 R 2 > 1. Comparison of Eqs. (4.16), (4.17) illustrates 
that tube length fluctuations cannot properly be sepa- 
rated from tube destruction. Being of the same order of 
magnitude for t w T 2 , they strongly interfere in the full 
theory 

We finally note that the use of T3 as fit parameter, 
necessary to fit the results of the primitive chain model 
to our data, illustrates that tube length fluctuations in- 
fluence the efficiency of tube destruction for all times [6] , 
even for chains of lengths N/N e w 100. 

Internal relaxation can be analyzed by considering a 
subchain of length M in the center of an infinitely long 
chain. We thus consider the coherent structure function 



S c (q,t;M, 00) = ^e^(*)-^(°))^ 
j,k=i 



(4.18) 



Within time t, segment j is displaced along the tube by 
£ s n steps, and since the tube has a random walk confor- 
mation, the average over the paths of segment j yields 



e iq(r 3 (t)-r fc (0)) 



-<?\j+l 3 n-k\ 



This expression is to be averaged over the distribution 
V\{n,t) of the number n of particles that diffused over 
bead j. We note that in the situation considered here 
Vi(n,t) is independent of j, since the process is transla- 
tionally invariant along an infinitely long chain. We thus 
find 



M 



S c (q,t;M,oo)= £ V^n,t) £ e"«" \^~ k \ . 

n=-oc j,fc=l 

(4.19) 

Now shifting j +l s n — > j, we can rewrite this expression 



as 



+00 



M 



S c (q,t;M,oo)= ^ Pi(n,t)l ^ e"« ^ 



n= — oo 
M 



M+l s n i s n 

f\3-k\ _ \ ^ p -t\3-k\ 



J2 e-rii-fcl-^e- 

3 = M+1 j=l 



(4.20) 



Since V\(n,t) is normalized to 1, the first term in Eq. 
(4.20) yields the static structure function of the subchain. 
All time dependence results from the motion of the end 
pieces of the subchain and thus is due to the length fluc- 
tuations and the shift of the subchain in the tube. Proper 
interior relaxation within the subchain has no effects on 
the coherent structure function. 

We want to stress that this does not imply that our 
theory neglects interior relaxation. The basic relation 
(4.19) only exploits the fact that the segment j has car- 
ried out a random walk of isn steps. It does not imply 
that this walk is along the precise original configuration 
of the chain. Rather the path will have changed by the 
motion of other segments, i.e., internal relaxation. What 
our result shows is that the effects of interior relaxation 
average out in the structure function as calculated within 
the reptation model. This also explains the observation 
in Subsect. 4.B that the discreteness of the individual 
steps is unimportant numerically It yields only a correc- 
tion to the second term in Eq. (4.20), which by itself is 
very small in the time range t < 10 3 where discreteness 
corrections might show up. Deviations among theory and 
data which we will find for short times and wave vectors 
of the order of the inverse segment size therefore must 
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originate from effects not taken into account in the rep- 
tation model. 

In this context some further remark on the model of 
Rouse motion in a coiled tube seems appropriate. Con- 
cerning the dynamics along the tube, our hopping model 
is essentially equivalent to Rouse motion. Specifically, 
beyond microscopic times the motion of individual seg- 
ments approaches Rouse type motion, as has been shown 
in Ref. [16], Sect. 3. The difference found for the struc- 
ture function results from the imposed initial conditions, 
only. It reflects the fact that the structure function mea- 
sures correlations between two segments j, k and thus 
is influenced by the initial configuration of the subchain 
connecting j with k. 

To summarize, we have seen that in the limit N — > oo 
with q 2 Rg ~ q 2 N kept fixed, the coherent structure func- 
tion of the total chain varies only on scale T3. In this 
limit, S c (q,t;N) is given by the primitive chain model. 
For shorter chains, tube length fluctuations come into 
play. They change the effective reptation time and lead 
to an initial decrease of S c on time scale T 2 . Both effects 
are important for chains up to lengths N/N e w 10 2 , at 
least. Proper internal relaxation, however, averages out. 

With regard to the data presented in Sect. 5, we want 
to close this subsection with a brief discussion of the scat- 
tering S c {q, t; M, N) from the central subchain of length 
M <C N in a chain of finite length N. For small times 
(and sufficiently large wave vectors), we expect to see 
small fluctuations of the length of the part of the tube 
occupied by the subchain. Since, however, these fluc- 
tuations do not lead out of the original tube, they are 
much less efficient in reducing the scattering function 
than the tube length fluctuations relevant for the total 
chain: on average, up to times of the order of the re- 
laxation time T 2 (M) of the subchain, its position in the 
tube has not changed seriously. A definite decrease of 
S c (q, t; M, N) sets in for times large enough for the sub- 
chain to leave its original part of the tube. For very large 
chain lengths N, this may occur for times t <C T 2 (N). 
S c then decreases like i -1 / 4 , this law reflecting the t 1 ^ 4 - 
law (Eq. (1.1)) for the motion of an internal segment. 
For T 2 (N) < t < T 3 (N), the motion of the subchain is 
driven by particles created by the chain ends which leads 
to a decrease like i -1 / 2 (cf. Eq. (1.1)). Finally, of course, 
this behavior is cut off by tube destruction reaching the 
position of the subchain. On the quantitative level this 
behavior can be extracted from our theory. In the dis- 
cussion of our data we in some figures will include curves 
which show the effect of relaxation of the (sub-) chain 
moving in an infinitely long tube, i.e., neglecting all tube 
renewal effects. 



model, can be solved only numerically. Furthermore, also 
the kernel and the inhomogencity in these equations do 
not take a simple form [3] . In particular, the inhomogenc- 
ity is given by some integral with a lengthy integrand. 
Thus, to present the results in a more easily tractable 
form some more or less heuristic analytical parameteri- 
zation might be useful. 

It should be clear that such a parameterization cannot 
be very simple. The coherent scattering function involves 
two time scales: T 2 and T 3 . Furthermore also wave vec- 
tor and chain length combine into two relevant variables: 
q 2 R 2 g and q 2 £ 2 y/N7N ~ q 2 R 2 ^jN e /N. The latter vari- 
able should govern the contribution of tube length fluc- 
tuations (cf. Eq. (4.15)). In searching for a reasonable 
parameterization we exploit the fact that for very long 
chains our results reduce to those of the primitive chain 
model. Furthermore, as shown in Sect. 3. A, the func- 
tional form of S c , given by this model, can be used to fit 
the long-time tail also for non-asymptotic chain lengths. 
This suggests to use an ansatz similar to de Gennes' form 
(3.9): 



SFi t {q,t;N) = (l-B)H(q,t;N) + BS DE \q 2 R 2 g ,^- 




The function Sde(Q, t) is given by Eqs. (3.2)-(3.5). Ana- 
lyzing the numerical results of our theory in a large range 
of wave vectors and chain lengths which exceeds that used 
in the simulations, we find that the effective reptation 
time weakly depends on q. We use the expression 



>,Flt 4 e \N P 
2.1 



2.5 



(9 + N/N e )0-5 9 + N/N e 



0.25 



(4 



(l + 0.3g 2 i? 2 ) 2 



The leading term gives the asymptotic law (4.12), the 
first correction is of order N e /N, in agreement with a 
result by Doi [6]. The prefactor of this correction dif- 
fers somewhat from Doi's result. However it must be 
noted that eqs. (4.21) - (4.25) are meant to reproduce 
the scattering function in a wide parameter range. Thus 
eq. (4.22) should not literally be compared to Doi's re- 
sult. 

We now turn to the coefficient B, which is taken as 



B = exp < 



q R g 



0.6 



2.8 



l + q 2 Rh 



D. Approximate analytical parameterization of our 
results for the total chain 

The integral equations which yield the normalized co- 
herent structure function S c (q, t; N) of the full reptation 



(4.23) 

This form is consistent with the discussion of the influ- 
ence of tube length fluctuations gi ven in t he previous sub- 
section. B tends to 1 for q 2 R 2 y / N e /N <C 1 and yields 

B»l- 0(yjN e /N) in the limit N -> 00, q fixed. 
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The initial decrease of the scattering function over a 
sizeable range of time is approximated by a stretched ex- 
ponential with an effective exponent which strongly de- 
creases with increasing q 2 R g - We take 



H(q,t,N) = cxp' 



T 2 



'-3,Fit 



3.2 



q 2 R% y t 

0.7 + q*R 2 g J f 2 
(4.24) 



xi = 0.1 + 



x 2 = 0.25 + 



0.9 



+ 0.2q 2 R 2 g 

0.75 
l + 0.1q 2 R 2 g 



(4.25) 



We recall Eq. (4.9): f 2 = (N+ I) 2 /it 2 . 

The parameterization (4.21)-(4.25) is fitted to the nu- 
merical evaluation of our theory in the range 0.1 < 
q 2 R 2 g < 50; 10 < N/N e < 340; < i < oo. In most 
of that range it reproduces our full results within devia- 
tions less than 0.015 in absolute value. It becomes worse 
for shorter chains (N/N e < 30) and larger wave vectors 
(q£o > 2) in the large time regime (t > O.IT3). In this 
region the large-time tail predicted by the full theory de- 
viates somewhat from the functional form given by the 
primitive chain model, and also small differences between 
the continuous and the discrete version of the model start 
to show up. For N/N e < 10 our theory breaks down, 
mainly due to the fact that tube length fluctuations be- 
come important even for t ~ T3. They thus affect the 
kernel V* {joto\m) of the integral equation (4.2), and this 
effect has not been considered in our theory. 

In the sequel we compare our Monte Carlo data to the 
results of the numerical evaluation of our full theory. It 
should be noted, however, that on the scale of the plots 
shown the difference between the results of the full theory 
and the effective parameterization (4.21)-(4.25) would be 
almost invisible. 



V. COMPARISON OF MONTE CARLO DATA TO 
FULL REPTATION THEORY 

The reptation model in full detail treats the motion 
along the tube axis, but it does not account for all relax- 
ation processes in a realistic chain. In particular, motions 
'perpendicular to the tube', like the creation and decay 
of larger side branches, will be present for the Monte 
Carlo chain but are not treated properly in the repta- 
tion model. Since larger side branches emerge by fusion 
of hairpins, they in the framework of our model would 
introduce some kind of interaction of the particles. Scat- 
tering vectors of magnitude {«1 certainly resolve such 
effects which result in some additional relaxation on a mi- 
croscopic time scale. In fitting the data to the theory we 
therefore allow for a phenomenological amplitude A < 1, 



multiplying the normalized structure function S c (q, t; N) 
calculated for pure reptation. We typically find values 
A w 0.9 for q = 1, A w 0.98 for q = 0.5, and Awl 
for all smaller q measured. Thus the deviations among 
theory and data which are found for small times, rapidly 
X2 \iecrease with decreasing q. 

} For a first illustration, we in Fig. 7 show our results 
lor N = 80, q = 0.5024 (q 2 R 2 g = 3.32). Data points 
result from runs extending to 10 5 or 10 s Monte Carlo 
steps. The full curve gives the result of our theory, with 
the amplitude adjusted to A = 0.985. Long dashes rep- 
resent a fit to the primitive chain model, with Boe and 
T 3 = t T 3 (MC) taken from Table 1. The short dashed 
curve results by neglecting tube destruction and is cal- 
culated as the structure function of a chain of length 
N in an infinitely long tube. Since the primitive chain 
model takes care only of tube destruction, the difference 
among the two dashed curves is due to tube length fluc- 
tuations. Thus Fig. 7 demonstrates that the tube length 
fluctuations strongly influence the structure function up 
to times t « T 2 . For much larger times tube destruction 
is the dominant process. 

These results are typical for all chain lengths and wave 
vectors. Figs. 8 and 9 show our results for N = 160 
and 320, respectively. Beyond some short time regime, 
we find excellent agreement between theory and data, 
with only the amplitude A — A(q, N) taken as adjustable 
parameter. Comparing Fig. 9 to Fig. 4, (note the dif- 
ference of times scales: log 10 i = \og 10 t^ MC ^ — 1.167), 
we again see the effect of tube length fluctuations: for 
N = 320, the primitive chain model can fit the data only 
for t > 2OT2. For smaller times, tube length fluctuations 
become relevant, their effect increasing with increasing 
q, as expected. For N = 160, q = 0.1 (q 2 R 2 g w 0.267), 
scattering cannot resolve the structure of the tube and 
the effect of tube length fluctuations is suppressed. The 
corresponding curve shown in Fig. 8 essentially coincides 
with the result of the primitive chain model, provided we 
adjust the effective time scale T 3 . 

To gain additional information on the chain dynamics, 
we also measured and calculated the structure function 
of the central piece of length M in a chain of length N. 
Fig. 10 shows our results for M = 80, N = 320, q = 0.5. 
With these parameter values it is directly comparable to 
Fig. 7, where M = N = 80, q = 0.5024. The central 
piece evidently is not influenced by tube length fluctua- 
tions. In contrast to Fig. 7, the decrease of S c (q, t; M, N), 
in Fig. 10 must be due to relaxation within the initial 
tube, up to times where tube destruction reaches the sub- 
chain measured. Indeed, in Fig. 10 the structure function 
calculated neglecting tube destruction (broken line) over 
a large time range coincides with the results of our full 
theory and simulations. The time Tr(M,N) it needs for 
tube destruction to reach the subchain, can be estimated 
from the relation 



i S n max (T R (M,N)) = 1(JV - M) 



(5.1) 
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This time is indicated in Fig. 10 and is close to the time 
where the dashed line starts to de viate fro m our full re- 
sult. Recall from Sect. 4. A that £sn max (t) gives the av- 
erage length of the end-piece of the tube that has been 
destroyed up to time t. 

As pointed out in Sect. 4.C, the decrease of 
S c (q, t; M, N) seen in Fig. 10 in the range t < T R (M, N) 
is not due to relaxation within the subchain M consid- 
ered, but results from the internal dynamics of the to- 
tal chain which leads to fluctuations of the position of 
the subchain in the original tube. This is illustrated in 
Fig. 11, which for given N — 320 and q = 1.0 com- 
pares the scattering from central subchains of lengths 
M = 40, 80, 160, 320. We note some peculiar be- 
havior: as function of time, the normalized scattering 
function of the total chain M = N = 320 initially de- 
creases faster and eventually crosses results for internal 
pieces. This shows the efficiency of tube length fluc- 
tuations in comparison to the internal dynamics of the 
chain. According to the discussion of Sect. 4.C, sub- 
chains M < N - yj \N e N w 300 (for N = 320) arc not 
influenced by tube length fluctuations, and for a large 
time range scattering from such subchains is determined 
by the diffusive motion of the total subchain in the ini- 
tial tube. Since the diffusion coefficient of this motion 
scales as 1/M, the diffusion of the longer subchains is 
very slow. Taking place in the initial tube, it is also less 
efficient in reducing S c than tube length fluctuations or 
the tube destructing shift of the total chain. For these 
reasons, the normalized coherent scattering for M = 160 
stays above the result for the total chain up to times 
t T 3 (N). For larger times the theory predicts that as 
for the shorter subchains the normalized scattering from 
subchain M = 160 falls below the scattering from the 
total chain, but the effect is too small to be visible in 
Fig. 11, and it clearly is not resolved by the accuracy 
of the data. We also note that tube destruction reaches 
the initial position of the subchain M — 160 only for 
pT^(160,320) w 10 5 - 8 . With decreasing M, diffusion in 
the tube becomes faster, and for the shortest subchain, 
M = 40, we see indications of a new regime: the slope of 
the curve seems to change near t = 10 5 4 . A doubly log- 
arithmic plot reveals that in the range 10 5 4 < t < 10 61 
the data approach a i _1 / 2 -law, and the upper bound of 
this interval is close to the time pTr(40,320) w 10 6 4 . 
These observations are consistent with the discussion of 
Sect. 4.C. 

Finally we consider the deviations between theory and 
data present for i < 10 3 and large wave vectors. A first 
important observation is that these deviations roughly 
are of the same size for interior pieces of the chain as 
for the total chain, see Fig. 11. Therefore their domi- 
nant source cannot be searched in our treatment of tube 
length fluctuations in the mean hopping rate approxima- 
tion. To show this in more detail, we in Fig. 12 have plot- 
ted the deviation between data and theory for wave vec- 
tor q = 1.0 and all chains and subchains measured. We 



note some end-effect: for the total chain N = 160 = M 
the deviation is about 20 % larger than for the central 
piece M — 160 in the chain N — 320. However, the 
bulk of the effect is independent of N > M, as illus- 
trated by the curves for M = 40 and M — 80, which 
each in fact combine data for N = 160 and N = 320. 
For these internal pieces, tube length fluctuations are ir- 
relevant, and in the time range shown, our evaluation 
of the reptation model is exact within the accuracy (0.5 
%) of our numerics. Consequently the initial deviations 
must be dominated by some relaxation in the interior of 
the subchain M which is not taken into account by the 
reptation model. We note that the amplitude A of the 
effect depends on M only weakly, but strongly increases 
with q, as is illustrated by Figs. 8 and 9. It roughly 
is proportional to q 2 . These observations are consistent 
with an interpretation as micro-structure effects, which 
clearly increase strongly with the resolving power of the 
experiment and are expected to show a weak dependence 
~ 1/M on (sub-)chain length. 

In our previous work on the motion of individual seg- 
ments, we found initial transients which also decayed on 
a time scale t < 10 3 . These could be traced back to 
the discreteness of both the hopping process and the un- 
derlying chain and could be largely explained by a fully 
discrete analysis of our model. As mentioned in Sect. 
4. A, it turns out that for S c (q,t, M, N) this discreteness 
numerically is irrelevant. In the light of the discussion 
of Sect. 4.C, this is quite understandable: the internal 
dynamics of the subchain averages out. We thus have 
to search for effects going beyond our implementation of 
the reptation model. One source of discrepancy immedi- 
ately comes into mind. A particle sitting on a bead is no 
exact representation of a hairpin. It, for instance, is not 
clear whether the bead sits in the tip or in the base of 
the hairpin. Simple estimates of such effects easily yield 
a correction amplitude 1 — A of the desired order of mag- 
nitude. However, the dynamics of single hairpins cannot 
explain a time scale of order 10 3 since it essentially re- 
laxes within one Monte Carlo step. Slower relaxation 
occurs for larger side branches. Already a structure cre- 
ated by fusion of two hairpins shows a relaxation time of 
the order of 10MC steps. Indeed, the fusion of hairpins 
destroys mobile units of the Monte Carlo chain. In the 
language of our model this induces some dynamical in- 
teraction of the particles. We therefore believe that the 
observed initial transients are due to the interference of 
reptation with the dynamics of larger side branches. 



VI. CONCLUSIONS 

In this paper, we presented results of simulations of 
the Evans-Edwards model, measuring the coherent struc- 
ture function of a rcptating chain in a large range of 
time and wave vector. This model contains the full rep- 
tational dynamics of a single chain, including internal 
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relaxation, tube length fluctuations, and tube renewal, 
but omits all the dynamics of the environment. In par- 
ticular, it ignores constraint release. We compared the 
data to two simplified versions of the reptation model as 
well as to our recent evaluation of the full theory. We 
found that the primitive chain model proposed by Doi 
and Edwards [4,2] yields a good fit to the data for times 
large compared to the Rouse time T 2 , provided we allow 
for a phenomenological amplitude factor and treat the 
reptation time T3 as a second adjustable parameter. De 
Gennes' approach [5] which adds a contribution meant to 
take one dimensional internal relaxation along the tube 
into account, cannot explain the data consistently. So 
the success of this approach in fitting neutron scattering 
data on polymer melts [21-23] indicates that the local 
contribution S^ e \q,t) to De Gennes' structure function 
phenomenologically accounts for some dynamics not con- 
tained in the reptation model. In view of our findings for 
the Evans-Edwards model we would search the origin of 
the observed initial decay in a coupling of reptation to 
motion perpendicular to the tube. Since in the experi- 
mental time range the initial decay accounts for most of 
the variation of the experimental curves [23] , we have not 
tried to fit these curves with our theory. 

Provided we allow for a single adjustable amplitude 
factor we, beyond the short time range influenced by mi- 
crostructure relaxation, find quantitative agreement be- 
tween our simulation data and our analytical evaluation 
[3] of full reptation theory. The theory in all details ex- 
plains not only the scattering function of the total chain 
but also the scattering from interior subchains. Except 
for the phenomenological amplitude, which notably dif- 
fers from 1 only for wave vectors of the order of the in- 
verse segment size, all parameters of the theory essen- 
tially are fixed by our previous analysis [16,17] of the 
motion of internal segments. Thus our version of the 
reptation model consistently explains both details of the 
internal motion and global features like the scattering 
function. 

Our analysis shows that the effects of internal relax- 
ation contained in the reptation model and seen, for in- 
stance, in the motion of internal segments, for the scat- 
tering functions average out. In the reptation model, 
the time dependence of the scattering functions is de- 
termined solely by the interplay of tube destruction and 
tube length fluctuations. The latter not only determine 
the scattering function up to times of order T2, but also 
for intermediate chain lengths lead to the decrease [6] of 
the reptation time T3 in comparison to the asymptotic 
law, clearly seen also in the scattering function. Numer- 
ically evaluating our theory, we found these effects to be 
well visible up to chain lengths of the order of 300 en- 
tanglement lengths. This is obvious from the numerical 
parameterization of our results, given in Section 4 D. 

Our theory deviates from the results of the simula- 
tion in the initial time region t < 10 3 . We stress that 
this is not due to the approximations inherent in our 
treatment of tube destruction. Deviations of the same 



order of magnitude occur for the total chain and for 
internal subchains, and the coherent structure function 
of internal subchains in the initial time region is not 
affected by our approximations. Thus the need to in- 
troduce some phenomenological parameter to take short 
time relaxation into account illustrates that even for 
the Evans-Edwards model, which is the most accurate 
computer-experimental implementation of the reptation 
model we can think of, some relaxation 'perpendicular' 
to the tube shows up. Clearly, for a melt, such effects 
must be much larger. However, our analysis suggests 
some way to suppress these contributions. The scatter- 
ing from a chain of length N = M should be compared 
to the scattering from the central piece of length M in 
a chain of length N 2M. Up to small end effects, the 
difference of the two coherent scattering functions should 
be due to reptation alone, and can be calculated from 
our theory. To observe an effect, the experiments must 
cover a time range extending at least up to the Rouse 
time, with chain length being of the order N/N e > 20. 
Unfortunately, for physical experiments, this seems to be 
out of reach. As shown by the work of Ref. [10] it might 
be feasible in simulations. 
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N (MC) 


q 


log 10 T, (MC) (g,iV) 

o 1 U o \i' / 


B D e(q, N) 

i_J 1_J V M 7 / 


80 


0.5024 


5.90 


0.91 


160 


0.10 


6.91 


0.997 




0.50 


6.87 


0.89 




1.00 


6.82 


0.76 


320 


0.25 


7.85 


0.965 




0.50 


7.82 


0.88 




1.00 


7.80 


0.78 



TABLE I. Parameter values for the fit of BDESDE{q 2 R 2 g ,t/Ts,) to the data. Reptation time 7^ MC ' is measured in Monte 
Carlo steps. 
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FIG. 1. A series of subsequent chain configurations, illustrating the microscopic dynamics of the Evans- Edwards model (in 
its two-dimensional version). The crosses represent impenetrable obstacles. 
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b) 

FIG. 2. Data for the coherent structure function (N = 320, q = 0.5) as function of log 10 t^ MC ^ in the short time regime. 
S c (g, t^ MC ^; N) is normalized to the exact static structure function So(q, N) (Eq. (2.6)). a) 10 sets of data, each averaged over 
10 3 independent runs. The lines serve to guide the eye. b) The same sets of data as in a), but with S c (q,0;N) subtracted. 



18 




-2 -1.5 -1 -0.5 0.5 



logio(t/T 3 ) 

FIG. 3. S c measured for q 2 R 2 g ~ 3.3 as function of t/Tz in the long time regime. Heavy dots: N — 320; small dots: N = 80. 
Full curve: primitive chain result (3.3). Ts(N = 320) has been adjusted to bring the long-time tail of the data for N = 320 on 
top of the theoretical curve. T 3 (80) = 73 (320) /4 3 . Some typical error bars (two standard deviations) are shown for N = 320. 
The statistical error rapidly decreases with decreasing t/T 3 . For tjTz < 10~ 0,5 it is smaller than the size of the dots. 
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FIG. 4. Data for S c (q, t, N) for chain length N = 320, wave vectors |q| = 0.25, 0.5, 1.0, as function of the Monte Carlo time; 



dots: medium time runs (tmi? = 10 8 ); circles: long time runs (t^ix 1 = 4.5 • 10 9 ). The curves give the results of the primitive 
chain model, fitted to the data as explained in the text. The arrow points to the equilibration time T^ 40 ^ (320), the heavy bar 
gives T3 (MC) (320) = 10 7 ' 82 . Some typical error bars (two standard deviations) are shown for q = 1.0, 0.25. 
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FIG. 5. Results for the reptation time defined in terms of the long time behavior of S c (dots) or in terms of the motion of the 
endsegment (circles). The asymptotic behavior T3 ~ TV 3 has been divided out. The lines correspond to power laws T3 ~ iV 3 ' 2 
(full line) or T 3 ~ N 3 5 (broken line). The short dashed line gives the asymptotic limit. 
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2 4 6 8 

FIG. 6. Fit of S dG (Eq. (3.9)) to data for N = 320, q = 1.0. Solid line: N e = 3.69, corresponding to B dG = 0.903; dashes: 
BdG = 0.82, ti (Eq. (3.8)) divided by 200. The arrow points to the equilibration time T^ ^ 1 . Data, small dots: short time 
runs (tmfjf' — 10 s ); heavy dots: medium time runs (tmfjf' — 10 8 ), circles: long time runs (t^l^ = 5 • 10 10 ) . 
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FIG. 7. Monte Carlo data for S c (q, t, N) compared to reptation theory (full line), as function of log 10 t. 
wave vector q — 0.5024 (q 2 R 2 g = 3.32). Long dashed line: best fit to the primitive chain model (Sde)- 
calculated neglecting tube destruction. The arrow points to Tb, the heavy slash gives T3. 
short time runs; heavy dots: averages over 50 medium time runs. 



Chain length N = 80; 
Short dashed line: S c 
Data, small dots: average over 10 4 
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logioi 

FIG. 8. Data compared to reptation theory. N = 160; q = 0.1 (q 2 R 2 g = 0.267); q = 0.5 (q 2 R 2 = 6.665); 
q = 1.0 {q 2 R 2 g = 26.61). Arrow: f 2 ; slash: f 3 . Small dots: short time runs; heavy dots: runs for intermediate time; 
circles: long time runs. The time regimes overlap. Some error bars (two standard deviations) are also given. 
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FIG. 9. Same as Fig. 8, but for N = 320; q = 0.25 (q 2 R 2 = 3.33); q = 0.5 (q 2 R 2 g = 13.33); q = 1.0 (q 2 R 2 g = 53.22). 
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FIG. 10. Coherent structure function of the central piece of length M = 80 in a chain of length N = 320, normalized to 
the static structure function of a chain of length M. Wave vector q — 0.5. Full line: full result of the theory; broken line: 
result neglecting tube destruction. The long slash through the curve at log 10 1 « 6.27 gives the time where tube destruction on 
average reaches the subchain M = 80, as derived from our theory. Other symbols as in Figs. 8 and 9. 
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FIG. 12. Deviation between theory and experiment: 8S = S^Monte Carlo) — A ■ S (theory), for q = 1.0. Connected data 
points: results for M = 40 (N = 160, 320), M = 80 (N = 160, 320), M = 160 (N = 320), (from below). For M = 40, 80 the 
results for the two different lengths N of the total chain fall right on top of each other. Disconnected points: M = N — 160. 
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